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the  box  causing  notable  bulging.  There  are  reports  of  a  sudden 
pressure  Increase  at  approximately  8000  ft/sec.  The  study  was  con¬ 
ducted  through  literature  search  and  simple  engineering  approximate 

The  study  concentrated  on  physical  phenomena  associated  with  ar 
impact  velocity  of  approximately  8000  ft/sec.  Heyda's  plug-shatter¬ 
ing  model  was  adapted  to  show  that  under  certain  impact  geometry  cor 
ditions  that  numerous  microsized  fragments  could  be  released  into 
the  closed  system,  at  or  above  the  impact  velocity.  At  approximate¬ 
ly  8000  ft/sec  these  microsized  fragments  encounter  stagnation  temp¬ 
eratures  slightly  greater  than  the  melting  point  of  aluminum  oxide, 
2300  K.  At  this  velocity,  the  stagnation  heating  begins  to  cause 
oxygen  in  the  air  to  dissociate  to  atomic  oxygen,  a  very  aggressive 
oxidizer.  These  conditions  are  compared  to  existing  experimental 
evidence  to  show  that  there  is  a  strong  likelihood  that  the  second¬ 
ary  fragments  would  ignite  and  burn  almost  instantly  releasing,  for 
a  multiple  impact  case,  a  significant  amount  of  energy  into  the 
system. 

The  Heyda  plug-shattering  model  was  extended  to  spherical  pro¬ 
jectiles.  Graphs  for  the  application  of  the  model  to  other  than 
alumlnum-alumir.um  systems  are  provided. 

It  is  recommended  that  further  experimental  confirmation  of  the 
predictability  of  the  synergistic  damape  phenomenon  be  performed 
before  direct  application  is  attempted. 
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Preface 


While  visiting  the  Air  Force  Armament  Laboratory  (AFATL) 
in  the  Spring  of  19?3»  the  results  of  an  experiment  were  pre¬ 
sented  in  an  overview  briefing.  The  first  picture  was  that 
of  an  aluminum  box  -like  structure  that  showed  damage  in  the 
form  of  multiple  punctures  in  both  the  front  and  back  sur¬ 
faces.  A  few  cracks  were  observed,  but  the  box  was  Intact. 

A  second  picture  of  a  similar  box  was  shown.  This  box  was 
more  severely  damaged;  there  were  large  cracks  and  evidence 
of  bulging  of  both  the  front  and  back  faces.  The  difference 
between  the  two  events  was  that  the  first  box  had  been  sub¬ 
jected  to  multiple  steel  fragments;  the  second  box  had  been 
subjected  to  multiple  aluminum  fragments.  Evidently  there 
was  a  mechanism  working  on  the  Al-Al  system  that  was  not 
present  in  the  steel-Al  system. 

Subsequent  conversations  with  AFATL  personnel  indicated 
that  there  appeared  to  be  a  "magic  number"  for  the  Al-Al 
system  -  8000  ft/sec  Impact  velocity  and  that  similar  behav¬ 
ior  had  been  observed  with  other  materials  impacting  A1 
structure.  The  cause  of  the  synergistic  damage  mechanism  was 
not  understood. 

The  purpose  of  this  study  was  then  to  present  a  reason¬ 
able  case  for  the  cause  of  the  synergistic  damage  by  focusing 
on  the  "magic  number"  -  8000  ft/sec.  I  must  express  my  thanks 
to  the  faculty  and  students  of  the  Graduate  Air  Weapons  class 
as  they  patiently  listened  to  dozens  of  hypotheses — and 
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explanations  why  that  hypothesis  wouldn’t  apply.  Especially 
my  thanks  to  my  advisory  committee,  Dr  Peter  J.  Torvlk,  MaJ 
Louis  Montulll  and  Capt  Wesley  Crow  who  bore  the  brunt  of 
these  discussions. 
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I,  Introduction 


Purpose  and  Scope 

The  purpose  of  this  study  was  to  hypothesize  the  most 
probable  cause  for  synergistic  damage  effects  observed  when 
multiple  aluminum  fragments  Impact  a  box-llke  aluminum  struc¬ 
ture  with  a  velocity  range  of  2-5  kra/sec.  The  study  was 
conducted  through  literature  search  and  simple  engineering 
approximations.  No  experimental  verification  was  contem¬ 
plated  since  the  Air  force  Armament  Laboratory  (AFATL)  would 
be  conducting  a  similar  study  based  on  a  comprehensive  test 
program.  Simple  engineering  design  data  criteria  were  to  be 
tentatively  established  for  the  Al-Al  system  and  some  exten¬ 
sion  made  for  other  materials. 

Background 

In  the  late  1940's  it  was  observed  that  when  high 
velocity  fragments  perforate  the  skin  and  Internal  structure 
of  an  aircraft,  flashes  or  semi-explosions  occur  which  may 
cause  considerable  bulging  and  tearing.  The  physical  nature 
of  the  phenomenon  was  attributed  to  burning  of  aluminum 
particles,  the  source  of  which  was  not  specified  (Ref  111). 

During  the  1960's  Investigations  of  hypervelocity 
meteorite  impacts,  several  investigators  reported  significant 
Impact  flash  when  aluminum  impacted  aluminum  plates,  espec¬ 
ially  at  velocities  greater  than  12,000  ft/sec  with  an  air 
or  oxygen  atmosphere  behind  the  plate.  The  Impact  flash  was 


1 


GAW/MC/74-10 


generally  attributed  to  the  burning  of  the  secondary 
particles  of  aluminum  (Ref  2,3). 

Backman  and  Stronge  reported  on  a  careful  Investiga¬ 
tion  of  the  Al-Al  Impact  system  but  were  unable  to  show  a 
significant  contribution  to  pressure  Inside  of  a  box  by 
combustion  of  A1  In  air.  They  show  a  general  correlation 
between  pressure  and  the  kinetic  energy  of  the  Impacting 
system  (Ref  4). 

More  recently,  work  sponsored  by  AFATL  has  shown  a 
significant  synergistic  damage  mechanism  when  multiple 
steel  fragments  impact  box-like  A1  structures.  The  velocity 
of  the  impacting  fragments  was  in  the  range  of  4000-8000 
ft/sec.1 

2 

Unpublished  information  from  AFATL  indicates  that  at 
an  Impact  velocity  of  approximately  8000  ft/sec,  for  some 
conditions,  a  significant  sudden  pressure  increase  can  be 
observed. 


^Extracted  from  unclassified  portion  of  a  classified  docu¬ 
ment..  For  specifics  contact  School  of  Engineering  ( ENE) 
Air  Force  Institute  of  Technology,  Wrlght-Patterson  AFB, 

Oh  45433. 

^Conversation  with  Dr.  McArdle,  AFATL/DLRD,  March  1973. 
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There  are  occasional  reports  of  sudden  pressure  Increase 
around  8000  ft/sec. 


roach 

At  least  two  observers  have  reported  a  sudden  increase 
of  pressure  or  Impulse  for  a  multiple  fragment  Impact 
velocity  of  8000  ft/sec.  By  understanding  the  phenomena 
associated  with  an  8000  ft/sec  Impact,  one  would  be  able 
to  deduce  a  physical  explanation  for  the  synergistic  damage 
effects  on  a  closed  system. 

First,  the  effect  of  multiple-fragment  Impact  was  con¬ 
sidered. 

Second,  to  understand  the  parameters  which  would  affect 
the  pressure  Increase  with  Impacting  kinetic  energy,  a  simple 
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drag  model  for  the  post-perforation  fragments  was  con¬ 
sidered. 

Third*  a  source  of  high-energy  secondary  fragments  was 
determined  by  a  plug-shattering  model  by  Heyda  (Ref  6) • 

Finally*  the  aerothermodynamlc  and  chemical  effects  on 
the  high  velocity  secondary  fragments  was  considered. 

Throughout  the  study*  factore  which  would  influence 
design  were  observed. 
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II.  Factors  Contributing  To  Synergistic  Damage 

This  section  will  summarize  the  significant  substudies 
which  lead  directly  to  an  explanation  of  the  damage  mechan¬ 
ism.  Many  other  factors  were  considered,  of  course,  but  are 
not  reported  either  because  It  was  believed  that  they  would 
lead  to  an  Insignificant  contribution  or  the  complex  Inter¬ 
actions  would  obscure  the  basic  phenomena. 

Multiple  Fragment  Impact 

The  first  exercise  was  to  construct  a  simple  computer 
model  which  would  provide  an  understanding  of  what  areal 
density  of  fragments  that  would  produce  significant  Inter¬ 
actions.  Visual  evidence  of  ruptured  plates  has  Indicated 
large  tears  and  evidence  of  tension  failure  In  addition  to  a 
number  of  perforations  when  the  synergistic  damage  was  ob¬ 
served.-^  For  this  reason  circumferential  stresses  were  con¬ 
sidered  more  important.  The  model  was  constructed  as  follows. 

Consider  several  holes  of  radius  a  which  are  punched 
out  of  an  Infinite  plate.  Stresses  of  magnitude  are 
generated  at  the  edge  of  each  hole  by  a  fragment  punching 
the  hole.  Assume  also  that  the  stresses  at  any  distance 
Y*  from  the  center  of  the  hole  is  given  by 


3see  Footnote  1. 
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where  07  Is  some  critical  circumferential  stress 
which  will  cause  failure  of  the  plate 
n  Is  a  dissipative  parameter.  For  this  model 
n  is  at  least  2. 

Now  assume  that  more  than  one  hole  is  punched  at  vir¬ 
tually  the  same  time.  The  stress  at  any  point  was  assumed 
to  the  summation  of  the  stress  emlnatlng  from  each  hole. 

The  results  of  the  analysis  Indicate  that  for  more  than  a 
few  hole  diameter  spacing,  the  Interaction  between  punctures 
Is  negligible. 

While  this  analysis  Is  extremely  simple,  It  is  also 
very  conservative.  Since  stresses  are  tensors  and  there¬ 
fore  havn  a  directional  property,  they  do  not  add  directly 
except  on  the  line  Joining  the  centers  of  each  hole. 
Mlklowltz'  very  much  more  elaborate  and  accurate  analysis  of 
the  stresses  Induced  by  a  punched  hole  Indicated  that  the 
more  rapidly  the  hole  Is  punched  the  faster  the  stresses 
dissipated  (Ref  5) •  Rinehart  and  Pearson  (Ref  6)  Indicate 
that  with  Impact  velocities  above  2000  ft/sec  the  shear 
stresses  at  the  edge  of  the  hole  are  negligibly  small  since 
the  shear  stresses  act  more  like  hydrodynamic  stresses. 

Therefore  one  can  reasonably  conclude  that  except  for 
very  closely  spaced  multiple  impacts  the  Increased  damage 
effects  are  probably  not  due  to  additive  or  synergistic 
effects  between  multiple  fragment  impacts.  Since  the 
evidence  Indicates  there  Is  a  much  greater  chance  of  observ¬ 
ing  the  sudden  pressure  bulld-up  for  multiple  fragment 


GAW/KC/74-10 


Impacts ;  there  must  be  some  other  phenomenon  which  more 
probably  occurs  because  cf  the  many  simultaneous  Impacts. 

It  will  be  shown  below  that  for  certain  fragment  sizes  and 
plate  thicknesses  one  could  get  a  planar  shock  In  the  plate 
which  will  produce  numerous,  high  velocity,  secondary  micro¬ 
size  particles. 

Energy  Release  Due  to  Particle  Ablation 

Consider  now  the  high  velocity  secondary  particles 
leaving  the  impacting  fragment-plate  system  Into  a  closed 
system,  the  boi.  Considerable  attention  will  be  given  to 
the  source,  size  and  velocities  of  these  particles  in  a 
later  section.  At  this  time  Interest  should  be  focused  on 
Backman's  observation  that  the  pressure  In  a  box  increases 
with  the  kinetic  energy  of  the  impacting  fragment  and  to 
the  frequently  observed  post-perforation  flash.  In  addi¬ 
tion  to  the  direct  mechanical  effects  of  the  Impact  causing 
the  pressure  rise,  these  particles  are  being  subjected  to 
aerodynamic  heating  which  may  cause  some  of  the  aluminum 
particles  to  ablate,  depositing  energy  In  the  form  of  heat 
In  the  closed  system. 

The  following  development  Is  after  Backman  (Ref  4t32), 

^  u 

An  energy  balance  on  the  1™  particle  gives 

*{;  s  .  Q  (2) 

where  is  the  heat  flux  incident  on  the  particle  per 

unit  area. 
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S  Is  the  cross  sectional  area  of  the  particle. 

Q  is  the  energy  being  carried  away  by  the 
melting  material  per  unit  time. 

The  Incident  heat  flu*  Is  equal  to  the  stagnation  pres¬ 
sure  times  the  current  velocity 


where  ^  Is  the  ambient  gas  density 

\y  Is  the  current  velocity  of  the  particle. 

For  a  particle  of  diameter  D  the  rate  of  energy  de¬ 


posited  Is 


where  C ^  Is  a  heat  exchange  coefficient. 

The  heat  being  carried  away  by  melting  Q  Is  given  by 


Q  •-  «[c,  (  u-  i.)  + 


where  fY\  Is  the  mass  of  material  being  carried  away 
from  the  fragment  per  unit  time 
C[  Is  the  specific  heat  of  the  material 
7J*  the  melting  temperature 
I  o  the  Initial  fragment  temperature. 

I'W  Is  the  heat  of  fusion. 

Equating  the  incident  heat  to  the  energy  being  carried 

* 

away  and  solving  for  M  ,  j 

TT  C •  P<  O'  V  (6) 

rw  =  _ &JJ - r^T 

6(cpT-v  hv  3 
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where  T  is  equal  to  (  T**i  -  1  «•  ). 


I  |^k 

Letting  w  sj^  ,  changing  variable  of  integration  to 


distance  traveled  X  rather  than  time,  we  have 


j__  z  c* ft  -  2  J 

8(CfT+H0 


(7) 


To  integrate  Eqn  (7)  an  expression  for  \/(x)  Is  required. 
If  the  only  drag  forces  are  acting  on  the  particle  in  the  x- 
direction,  Newton's  second  law  gives 


m,5‘-£-gV'  ^ 


(8) 


where  C6  is  the  drag  coefficient 

Mp  is  the  mass  of  the  fragment 
>(  is  the  second  time  derivative  of  distance. 
Again,  changing  variables  and  integrating 

V(x)»V.«rl 


8  Hr 


(9) 


Substituting  (9)  into  (7), 


T  G  ft  t>1 14 


<U  =  "  exp 


8(c,r  +  i>») 


T  £>*0)  X 


(10) 


Integrating,  the  mass  AM  which  is  transferred  from  the 
rapidly  moving,  and  decelerating  particle,  to  the  closed 
system  in  a  distance  AS  is  given  by 


C.  V.  mr 

Am  =.  *  - 


C„(CP  »  +  0 


“  exfr 


3fL  Co  D  A5 


X 


(ID 


The  value  of  is  on  the  order  of  unity  for  conditions  of 

significant  mass  loss  due  to  ablation  (Ref  7*224). 
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Substituting  into  Eqn  (11)  the  equivalent  value  for 
the  particle  mass,  and  area  in  terms  of  particle  size  and 
density,  we  have  after  correcting  for  the  assumed  spherical 
fragment  shape. 


Ajn 

/yip 


c* 

3Ct(Crni 0 


(12) 


where  is  the  density  of  the  fragment. 

The  left  side  of  Equation  12  represents  the  fraction  of 
the  energy  transferred  from  the  kinetic  energy  of  the  particl 
to  the  surroundings.  This  number  is  in  general  a  very  small 
value  unless  D  is  very  small,  is  very  small  or  there  are 
a  large  number  of  particles  contributing  a  small  fraction  of 
their  energy  to  the  system.  Figure  4  represents  equation 
12  plotted  for  10,800  ft/sec.  The  calculations  were  termi¬ 
nated  at  Vi  *  4000- ft/sec  because  that  is  the  minimum 
velocity  at  which  the  stagnation  temperatures  are  great 
enough  to  melt  aluminum.  Notice  that  with  particle  sizes  in 
the  10-100  micron  size  the  contribution  is  several  orders 
higher  than  for  the  somewhat  larger  fragments. 

Several  investigators,  notably  Eackman  and  Stronge 
(Ref  4:9)  made  determined  efforts  to  collect  the  fragments 
after  moderately  high  velocity  impact  events.  For  1/8*' 
thick  plates  and  impact  velocities  of  3.16  and  3*7  km/sec, 
Backman  reports  that  only  44#  and  55#  of  the  total  mass  of 
the  projectile  and  target  were  recovered.  Most  of  the  frag¬ 
ments  recovered  were  fairly  large,  greater  than  ~  250 . 
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For  the  higher  velocity  impacts  a  considerably  higher  per¬ 
centage  of  the  recovered  fragments  were  quite  small,  the 
smallest  sieve  size  being  #200  (an  opening  of  74  ). 

From  these  data  and  equation  12  there  exists  the  sug¬ 
gestion  that  there  are  a  considerable  number  of  very  small 
fragments  moving  at  very  high  velocity  exiting  the  fragment 
hole  into  the  box.  Equation  12  also  suggests  an  analytic 
confirmation  of  Beckman' s  observation  that  the  pressure 
build-up  in  the  box  varies  with  the  kinetic  energy  of  the 
impacting  system. 

What  then  could  be  the  mechanism  for  the  sometimes  ob¬ 
served  sudden  pressure  increase  at  approximately  8000  ft/sec? 
It  is  the  hypothesis  of  this  study  that  under  certain  condi¬ 
tions  that  very  tiny  highly  energetic  particles  are  heated 
aerodynamically  to  the  ignition  temperature  and  that  extreme¬ 
ly  rapid  burning  occurs.  The  first  consideration  then  is  to 
show  a  way  of  predicting  under  what  conditions  these  second¬ 
ary  particles  occur.  It  will  then  be  shown  that  the  stagna¬ 
tion  temperature  at  approximately  the  critical  velocity  of 
8000  ft/sec  in  air  is  such  that  there  is  significant  disso¬ 
ciation  of  molecular  oxygen  to  atomic  oxygen.  Also  at  that 
velocity  the  stagnation  temperature  is  at  or  above  the  melt¬ 
ing  temperature  of  aluminum  oxide  (AI2O3)  allowing  rapid  re¬ 
action  of  the  energetic  oxidizer  and  aluminum. 

Heyda  Plug-Shattering  Model 

A  likely  source  of  energetic  microsized  particles  has 


12 


GAW/MC/7^-10 


been  hypothesized  by  Heyda  (Ref  9*1^2).  The  model  predicts 
that  under  certain  conditions  the  plug  driven  out  of  a 
plate,  subjected  to  high  velocity  Impact  by  a  flat-faced 
cylindrical  projectile,  will  shatter  at  or  near  the  center. 
For  the  few  cases  where  verification  was  attempted,  the  Im¬ 
pact  velocity  of  Incipient  plug-shattering  was  accurately 
predicted.  The  basic  development  of  the  model  Is  given  with 
reference  to  figure  3» 


Fig  3.  Geometry  and  Nomenclature  for  Development  of  Heyda 
Plug-Shattering  Model  (Ref  9*1^3)* 

At  time  ±T  the  projectile  plate  Interface  has  moved  with 
velocity  U  ,  the  particle  velocity  In  the  plate,  a  short 
distance  Into  the  plate  with  thickness T  .  Figure  3  repre¬ 
sents  the  conditions  at  tj  ,  the  Instant  that  the  planar  shock 
front  traveling  at  velocity  D  Just  reaches  the  rear  surface 
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of  the  plate.  The  compressed  region  behind  the  shock  wave 
Is  being  attenuated  by  a  relief  wave  eminating  from  tie 
corner  of  the  projectile  at  C.  This  relief  wave  move*  at 
the  hydrodynamic  sound  velocity  C„  for  the  conditions  behind 
the  shock.  Observe  that  the  center  of  the  rarefaction  wave 
is  moving  at  the  particle  velocity  behind  the  shock  U  . 

At  the  instant  the  shock  reaches  the  rear  surface,  the  relief 
wave  has  started  to  relieve  the  intense  pressure  at  B.  De¬ 
noting  the  radius  of  the  unattenuated  shock  A£  as  X,  then  at 
time  tT  the  following  hold. 

0  -  w  =  Cg  a 

(13) 

R  -  X  =  (c*  cos  *  )  ^ 


T  =  Dir 


(14) 

(15) 


combining  these  equations 


and  after  considerable  algebraic  manipulation 


X  = 


(16) 


(17) 


where 

v/c;  -  (  l  -  “X  (1 

Note  that  parameter  S  is  a  dimensionless  function  only  of 
the  material  properties  and  the  huponiot  conditions  which 
are  specified  by  the  initial  velocity. 
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Referring  to  Figure  3  and  equation  17  observe  that  if, 

S  >  ft  =7*  X  >0  (19) 

and  the  shock  wave  will  not  be  attenuated  on  the  arts  of 
symmetry  before  reaching  the  rear  surface  of  the  plate,  A 
reflection  of  the  compression  wave  will  occur  at  the  air- 
metal  interface  and  the  familiar  spallation  phenomenon  may 
occur  if  the  reflected  tensile  wave  is  of  sufficient  strength 
and  shape. 


Fig  4.  Plot  of  Dimensionless  Variable  S  Versus  the  Ratio  of 
Material  Density  Before  Impact  to  Density  After  Impact 
(Ref  9*144).  Additional  Plots  for  a  Wide  Variety  of 
Materials  are  Shown  in  Appendix  B. 
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Now,  consider  again  Figure  3  and  equation  17. 


If 

S  <  "jjr  =>  x  <o 


(20) 


which  implies  that  the  shock  wave  will  be  attenuated  on  the 
axis  of  symmetry  before  reaching  the  rear  free  surface. 

"Such  an  attenuation  point  corresponds  physically  to  the 
rapid  deposition  of  energy  on  the  ails  of  penetration  at  the 
point  which  Is  T-RS  units  from  the  rear  surface.  This 
rapid  energy  deposition  approiimates  eiplosive  behavior  and 
results  in  plug- shattering. "  (Ref  8*143) 

From  equation  19  and  20  one  concludes  that  there  exists 
a  critical  velocity  V4  >  at  or  above  which,  plug-shattering 
is  expected  to  occur.  The  critical  case  occurs  when 

5  *  "R  (21) 
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Fig  5»  Case  Where  S  * ^  .  Energy  is  Deposited  at  Plug 
Centerline  af,  an  Explosive  Rate. 


To  calculate  the  critical  velocity  Va  at  which  incipient 
plug-shattering  occurs  begin  with  the  Hugoniot  continuity 
equation  across  a  shock  wave 


p.D  =  p(0-u) 


(22) 


where  p.  is  the  ambient  material  density,  p  is  the  density 
of  the  material  Immediately  behind  the  shock,  [)  is  the  shock 
front  velocity  and  U  is  the  particle  velocity  in  the  material 
behind  the  shock. 

Combining  equation  22  with  the  well  known  eraperical 
relation 


D  =  A+ 


(23) 
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where  A  and  B  are  emperical  constants  the  values  of  which 


are  summarized  in  AFWL-TB-69-38  (Ref  9) ,  one  gets 


A 


for  the  particle  velocity  in  the  plate. 


(24) 


In  the  projectile,  equation  23  becomes 

5  =  a  +  b(^«-w) 


(25) 


where  the  superscript  bar  implies  properties  of  the  projec¬ 
tile. 

Finally,  after  considerable  manipulation,  it  can  be 
shown  (Ref  8:145)  that  the  critical  velocity  for  incipient 
plug-shattering  is  given  by 

~  A  (26) 

where  is  the  density  of  the  projectile  prior  to  im¬ 

pact* 

Experimental  data  contained  in  Woodall  (Ref  8:146)  shows 
for  a  steel  cylinder  of  radius  0.391  in  impacting  a  steel 
plate  of  thickness  0.635  In  at  a  velocity  of  .777  km/sec 
(2550  ft/sec)  showed  signs  of  shattering  type  failure  at  the 
plug's  axis  of  symmetry.  A  calculation  of  the  expect’ed 
critical  velocity,  the  velocity  at  which  incipient  plug-shat¬ 
tering  would  occur,  indicates  Vo  should  be  2619  ft/sec.  This 
suggests  excellent  agreement  between  experiment  and  a  simple 
model.  The  details  of  the  calculation  are  shown  in  Appendix 
A. 
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Extrapolation  of  the  Heyda  Plug-Shattering  Model  to  Spherical 
Projectiles 

Above  we  have  seer,  that  the  Heyda  plug-shattering  model 
predicted  fairly  well  when  Incipient  shattering  will  occur 
for  a  flat-faced  cylindrical  projectile  Impacting  a  plane 
plate.  Since  It  Is  stated  by  Swift  (Ref  10*518)  that  little 
is  known  quantitatively  about  the  material  shattering  pro¬ 
cesses,  it  may  be  appropriate  to  make  a  heroic  extension  of 
Heyda* s  model  to  other  shapes. 

Heyda 's  derivation  assumes  a  planar  shock  Is  generated 
by  the  plane  face  of  the  cylinder  which  creates  a  region  of 
very  high  pressure  in  a  region  behind  the  shock  which  Is 
then  relieved  by  a  rarefaction  wave  converging  on  the  axis 
of  symmetry  with  explosive  results. 

Consider  a  sphere  impacting  a  plate.  At  the  point  of 
impact  at  the  instant  of  impact  there  will  be  some  deforma¬ 
tion  of  the  projectile  and  plate  which  would  approximate  a 
planar  impact  of  a  smaller  cylinder.  Note  that  Heyda* s 
model  depends  only  on  material  properties,  Impact  velocity 
and  geometry  of  the  Impact,  there  is  no  restrictions  on  the 
mass  (or  length)  of  the  projectile. 

To  make  such  an  extension  of  Heyda *s  model  one  would, 
of  course,  prefer  to  predict  on  sound  theoretical  grounds 
which  parameters  would  or  could  equate  different  shapes.  In 
this  extension  for  example,  an  equivalent  radius  R*  Is  pro¬ 
posed  which  equates  a  spherical  contact  surface  to  a  plane 
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circular  surface.  It  appears  to  be  a  formidable,  If  not 
Impossible,  task  to  find  such  a  theoretical  relation.  A 
second  choice  would  be  to  use  one  of  the  very  complex  (and 
time  consuming)  computer  hydrodynamic  codes  such  as  CRAM, 

i 

PICKWICK,  VISTA  OR  KELP  (Ref  ll»l63).  An  objective  of  this 
study  was  to  minimize  complex  calculations  for  the  first 
order  analysis  so  the  computer  codes  were  not  extended  to 
include  this  problem.  The  final  alternative  was  to  look  at 
the  very  limited  experimental  data  which  reported  informa¬ 
tion  from  which  one  might  attempt  an  extension.  The  best 
source  of  these  data  are  reported  in  Hopkins,  Lee  and  Swift 
(Ref  12). 


Fig  6.  Geometry  of  Spherical  Projectile  Impacting  a  Plate 
for  Extention  of  Heyda  Plug-Shattering  Model. 
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IMPACT  VELOCITY,  km/scc 


Fig  7.  Ballistic  Limit  Curves;  Is  Thickness  of  Witness 
Plate  (Ref  12) 

Before  describing  the  method  for  estimating  R*,  the 
equivalent  radius  for  a  sphere  which  allows  one  to  use 
Heyda's  model,  an  explanation  of  the  source  of  empirical 
data  Is  required.  Hopkins,  Lee  and  Swift  were  attempting 
to  learn  more  of  the  role  of  shock  heating  in  determining 
the  performance  of  two  element  hypervelocity  impact  shields. 

A  series  of  tests  were  run  where  an  aluminum  projectile  of 
3.17  mm  diameter  was  propelled  through  an  aluminum  bumper 
plate  0.787  mm  thick.  The  thickness  of  the  hull  plate  con¬ 
sisting  of  606I-T6  aluminum  was  varied  until  a  ballistic 
limit  was  found  for  a  given  impact  velocity.  Figure  7  shows 
the  thickness  of  the  hull  plate  at  the  ballistic  limit  for 
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a  wide  range  of  Impact  velocities.  At  the  point  Indicated 
as  "Incipient  shattering"  some  part  of  the  system  Is  pre¬ 
sumed  to  shatter  Into  smaller  pieces,  assumed  to  be  the 
plug  In  accordance  with  the  Heyda  model. 

Using  3»0  km/sec  taken  from  Figure  7  equations  24  and 
22  were  solved  In  reverse  for  —  .  Using  figure  4,  or 

charts  shown  in  Appendix  Bc  a  corresponding  S  was  found. 

Since  at  the  critical  velocity  5=  T/r  and  1  was  fixed, 
an  equivalent  R*  was  found.  The  ratio  was  found  to  be 

between  .25  and  .30  for  all  experimental  evidence  of  incipient 
plug-shattering  found.  The  details  of  one  calculation  are 
shown  In  Appendix  A. 

In  summary,  an  argument  has  been  presented  which  Implies 
that  a  source  of  energy  which  causes  a  pressure  increase  In 
a  closed  box  Is  dependent  on  numerous,  very  tiny,  rapidly 
moving  secondary  fragments.  One  source  of  these  fragments 
under  certain  circumstances  of  velocity,  plate  thickness  and 
Impact  particle  size  and  density  would  be  the  explosive  re¬ 
lease  of  energy  in  the  center  of  the  plug  as  predicted  by 
the  Heyda  model.  Presumably  these  fragments  would  have  a 
velocity  greater  than  that  which  would  be  expected  from  a 
like-material  impact,  with  spallation.  For  a  spallation 
event  the  secondary  fragments  have  an  Initial  velocity  of 
twice  the  plate  particle  velocity  which  Is,  for  llke-material 
Impact,  one-half  the  Initial  Impact  velocity.  Or  more  con¬ 
cisely,  for  like-materlal  spallation  event  the  secondary 
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fragments  will  have  the  same  velocity  as  the  Impacting 
projectile. 

It  is  also  true  that  there  will  be  other  small  frag¬ 
ments  of  somewhat  lower  velocity  than  either  the  plug,  the 
impacting  projectile,  the  traditional  spallation  fragments 
or  the  shattered  plug  fragments  (If  the  conditions  exist). 
These  small  fragments  at  a  lower  velocity  are  associated 
with  the  buildup  of  a  combination  of  tensile  and  shear 
stresses  in  the  plate  induced  by  the  wedging  interaction 
between  the  projectile  and  target  as  the  projectile  passes 
through  (Hef  8:377), 

There  is,  therefore,  a  wide  distribution  of  velocities 
for  the  secondary  fragments  possible.  It  would  be  difficult 
if  not  Impossible  to  specify  an  exact  distribution  for  all 
cases.  For  the  purpose  of  further  analysis  microsized 
fragments  will  be  assumed  to  come  from  the  plug-shattering 
event.  Further,  the  assumption  will  be  made  that  at  least 
some  of  these  fragments  will  have,  as  a  minimum,  the  velo¬ 
city  of  a  spallation  fragment, 8000  ft/sec.  Complex  com¬ 
puter  code?  may  be  modified  to  give  a  more  accurate  estimate. 

The  next  step  is  to  consider  what  effect  these  frag¬ 
ments  may  have  on  the  pressure  build-up  in  a  box. 

Aero thermodynamic  Effects  on  Secondary  Fragments 

Having  considered  a  source  of  high  velocity  secondary 
fragments  and  focusing  attention  on  the  very  small  fragments, 
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in  the  tens  of  micron  size  and  a  velocity  of  8000  ft/sec, 
the  next  step  is  to  consider  the  aerothermodynamlc  effects 
on  these  fragments  and  their  immediate  surroundings.  It  Is 
the  hypothesis  of  this  study  that  these  fragments  will  burn 
almost  immediately  at  or  above  the  critical  velocity  being 
considered.  In  slmplist  terms  the  requirements  for  burn¬ 
ing  are*  (1)  a  fuel,  a  material  that  will  react  rapidly  with 
oxygen  (or  other  element)  exothermally;  (2),  an  oxidizer, 
oxygen,  in  the  form  of  atomic  oxygen,  0,  not  02  or  the  oxygen 
ion  0~;  and,  (3)  a  source  of  heat  which  will  raise  the 
temperature  of  the  system  to  the  ignition  temperature. 

Fragment  Heating.  Before  making  claims  as  to  the 
amount  of  heating  of  the  air  it  is  necessary  to  know  if  the 
flow  field  can  be  treated  as  a  continuum,  in  which  the  com¬ 
monly  used  hypersonic  analysis  have  been  performed.  The 
usually  accepted  criteria  for  continuum  flow  is  the  Knudsen 
number,  ,  which  is  the  ratio  of  the  molecular  mean  free 
path,  ,  to  some  characteristic  dimension  of  the  flow 
field  L  .  If  K*  is  very  much  less  than  unity  then  the  flow 
is  a  continuum  (Ref  13s 210).  At  sea  level  for  a  50  micron 
fragment , 

X.  30X10  ^  =  C,X/o*V  <<  [ 

L  So  X  «e  ^  (27) 

Therefore,  the  classical  continuum  analysis  will  be  used  for 
a  real  gas. 

The  first  consideration  is  the  temperature  of  the  air 
surrounding  the  fragments.  Figure  8  shows  the  temperature 
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Immediately  behind  a  hypersonic  shook  wave  of  a  cylindrical 
fragment  moving  In  air  at  sea  level  and  35.000  ft.  Note 
that  at  8,000  ft/sec  the  temperature  Is  approximately  2790  K 
If  the  initial  temperature  mas  300  K.  The  boiling  tempera¬ 
ture  of  A1  IS  2350  K  (Ref  15). 


Figure  9  represents  the  temperature  at  the  stagnation 
point  or.  a  cylindrical  fragment.  At  3000  ft/sec  the  tempera¬ 
ture  is  approximately  2730  K.  It  Is  well  known  that  the 
Ignition  temperature  for  aluminum  Is  approximately  2300  K 
(Hef  16)  which  is  below  the  ambient  temperature  of  a  particle 

O  with  a  velocity  of  8000  ft/sec. 
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Fig  9*  Stagnation  Temperature  for  a  Circular  Cylinder 
(Ref  15:163) 


One  would  then  be  concerned  why  all  aluminum  projec¬ 
tiles  do  not  burn  up  entirely  when  traveling  at  velocities 
in  the  neighborhood  of  8000  ft/sec.  Fay  &  Riddell  (Ref 
13:183)  point  out  that  the  heat  transfer  rate  to  the  body 
is  a  function  of  the  air  density  p  ,  and  the  radius  of  the 
entry  face  of  the  projectile  rb  .  Figure  10  represents  the 
relationship  between  these  parameters  and  the  free  stream 
velocity.  T[r-  is  the  temperature  of  the  body.  Note  parti¬ 
cularly  that  for  very  tiny  fragments,  micron  size,  that  the 
heat  transfer  rate  to  the  body  is  extremely  high,  In  fact 
the  smaller  the  fragment  the  faster  it  will  heat  to  the 
ignition  temperature  and  supply  the  necessary  heat  required 
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for  any  energy  threshold  which  might  exist  which  would  delay 
ignition.  The  key  point  on  Fig  10,  however,  is  that  at 
approximately  8000  ft/sec  for  a  very  hot  body,  the  heat 
transfer  rate  achieves  a  positive  value.  Ignition  will 
occur  only  if  enough  energy  is  added  to  the  system  to  heat 
a  slab  of  aluminum  on  the  fragment  as  thick  as  a  steadily 
propagating  adiabatic  laminar  flame  to  the  ignition  tempera¬ 
ture  (Ref  16:188) . 


Fig  10.  Stagnation  Point  Heat  Transfer  in  Laminar  Equili¬ 
brium  Dissociated  Flow  (Ref  13:184) 

Effects  on  Ambient  Atmosphere  Surrounding  Fragments. 

At  8000  ft/sec  the  secondary  fragments  are  exposed  to  temp¬ 
eratures  of  approximately  2700  K  or  400  K  greater  than  the 
ignition  temperature  of  aluminum.  This  section  will  address 
the  concentration  of  high  temperature  oxidizer  surrounding 
the  fragments. 
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Pig  12.  Stagnation  Pressure  at  Leading  Edge  of  a  Right 
Cylinder  at  Hypersonic  Velocities  (Ref  15:178) 


Another  consideration  is  the  chemical’ effects  on  air  at 


these  very  high  velocities.  Figure  13  suggests  a  very  im¬ 


portant  consideration,  coincidental  with  8000  ft/sec  flow. 


At  that  velocity  dissociation  of  molecular  oxygen  into  atomic 


oxygen  becomes  significant. 
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Fig  13.  Chemical  Effects  of  Air  at  Various  Velocities  and 
Altitudes  (Ref  1?:21) 


Another  representation  of  similar  data,  taken  from  Cox 
and  Crabtree  (Ref  13il6l)  relates  the  compressibility  factor, 
Z  which  is  defined  by 


Z* 


(28) 


where  !  is  the  gas  pressure 

p  is  the  gas  density 

$  is  the  universal  gas  constant 
T  is  the  absolute  temperature  of  the  gas 
yr\  is  the  molecular  weight  of  the  gas 
in  consistent  units.  The  parameter  Z  can  be  represented  by 
l  +#  ,  where  is  the  fraction  of  initial  moles  of  gas 
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which  has  dissociated  (Ref  3*108)  and  therefore  represents 
a  deviation  from  Ideal  gas  law.  The  greater  the  value  of 
the  more  the  gas  deviates  from  an  ideal  gas.  From  Fig 
14  note  that  Just  above  a  temperature  of  2500  K  the  dissocia¬ 
tion  of  molecular  oxygen  begins  to  change  the  behavior  of  a 
gas.  This  deviation  from  ideal  gas  behavior  has  been  Implicit 
in  all  arguments  presented  in  this  section. 


rtf 


TEMPERATURE  ,Kx  I0'a 


Fig  14.  Compressibility  Factor  for  Air  (Ref  18 a  123 ) 


Combustion  and  I gnl tlon  of  Aluminum  Fragments.  In  the 
past  two  subsections  it  was  shown  that  if  fine  particles  of 
aluminum  are  traveling  through  air  at  approximately  8000  ft/ 
sec  that  the  aero thermodynamic  heating  will,  coincidentiy , 
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cause  temperatures  surrounding  the  particle  to  rise  to  the 
Ignition  temperature,  and  that  the  smaller  the  fragment  the 
greater  the  heat  transfer  to  the  fragment.  Additionally,  at 
that  critical  velocity  air  begins  to  dissociate  to  atomic 
oxygen,  an  extremely  powerful  agent  for  supporting  burning 
of  aluminum,  or  any  other  combustible  material.  The  next 
subject  will  address  briefly  some  of  the  phenomena  necessary 
for  an  aluminum  particle  to  burn  almost  Instantaneously,  In 
the  environment  that  exists  around  the  fragment  while  moving 
at  approximately  8000  ft/sec. 

Friedman  and  Macek  (Ref  18»15)  in  their  theoretical 
and  experimental  Investigation  of  the  burning  of  aluminum 
particles  conclude  that  the  melting  point  of  alumina 
(A^C^)*  2300  ^  20  K,  Is  very  Important  for  ignition  of 
aluminum.  At  that  temperature  the  alumina  melts  which 
allows  the  diffusion  of  the  oxygen  to  the  aluminum.  They 
point  out  that  theory  predicts  that  the  Ignition  should 
depend  strongly  on  the  concentration  of  oxygen,  but  as  Is 
shown  by  Fig  15»  the  ignition  limit  is  relatively  insensi¬ 
tive  to  the  concentration  of  oxygen. 
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Fig  15.  Experimental  Ignition  and  Fragmentation  Limits 
(33  M  Diameter  Aluminum)  (Ref  IB) 

*  Minus  Concentration  Implies  Deficiency  From  Stoichio¬ 
metric  Concentration. 


The  key  to  combustion  and  combustion  rates  at  elevated 
temperature  then  seems  to  be  that  a  reaction  will  begin 
where  the  oxide  coating  surrounding  the  at  least  partially 
molten  fragment  can  be  penetrated  by  oxygen  or  as  in  the  case 
of  the  high  velocity  environment,  atomic  oxygen.  Almost 
instantly,  the  temperature  rises  to  greater  than  2350  K,  the 
boiling  point  of  aluminum,  and  the  reaction  takes  place  in  a 
gaseous  system.  For  observed  cases  an  increase  in  pressure 
is  recorded  approximately  50  yu  see  after  impact  (Ref  4i46-54) 
so  the  ignition  time  is  of  utmost  importance,  if  burning  is 
to  be  considered  an  important  contribution  to  the  sudden 
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pressure  increase.  Delayed  ignition  could  also  contribute 
significantly  by  increasing  the  duration  of  the  unaugoented 
pressure  pulse  associated  only  with  a  high  velocity  impact. 

Belyaev  (Ref  19023)  measured  the  ignition  time  and 
burning  time  of  aluminum  particles  in  "jh  pressure  oxidiz¬ 
ing  atmosphere.  The  pressure  range  a  .red  was  10-100 
atmospheres  which  was  shown  in  Fig  11  and  12  to  be  the  applic¬ 
able  case  for  8000  ft/sec  flow.  Their  results  show  that  for 
a  ?0 /j.  diameter  fragment  the  ignition  time  is  essentially 
zero  for  temperatures  above  2500  K.  Fig  16  illustrates  their 
data. 


Fig  16.  Ignition  Time  for  70/*.  A1  Particles  in  High 
Temperature  Air  (Ref  19024) 


s 
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They  also  show  that  the  burning  time  Is  a  strong 

function  of  the  particle  diameter  and  the  strength  of 

r 

oxidizer  which  they  characterize  by  the  parameter  Q .  . 

They  did  not  consider  the  effect  of  atomic  oxygen  In  the  re¬ 
action.  represents  the  activity  of  the  oxidizer  and  the 

concentration,  a  higher  value  Implies  a  more  energetic  re¬ 
action.  An  assumption  must  be  made  that  a  high  concentra¬ 
tion  of  atomic  0  would  be  a  more  energetic  atmosphere  than 
an  atmosphere  where  the  oxygen  carriers  are  CO,  H20,  or  C02 
as  was  the  case  In  which  their  experiments  were  carried  out. 

The  amount  of  energy  released  by  burning  fragments  must, 
therefore,  depend  on  the  burning  time.  Figure  17  shows 
Belyaev's  data  for  burning  time  of  aluminum. 
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Fig  18.  Burning  Time  of  A1  Particles  as  Function  of 

Activity.  a.£  Represents  Concentration  and  Availability 
of  Free  Oxygen,  Increasing  a.*  Implies  Greater  Activity. 
(Ref  20:323)  *  J 


From  Fig  17  and  19  it  is  clear  that  the  rate  of  energy 
release  is  a  strong  function  of  the  oxidizing  capability 
of  the  atmosphere  and  the  size  of  the  fragments.  As  was 
mentioned  above,  at  approximately  8000  ft/sec  the  tempera¬ 
tures  surrounding  the  fragment  are  approximately  2500  K  and 
that  air  begins  to  deviate  from  ideal  gas  behavior  because 
of  dissociation  of  0£  to  atomic  oxygen, 

A  parameter  which  will  Indicate  the  amount  of  0  in  the 
surroundings  is  the  equilibrium  constant  Kp  .  \(p  ,  based 

on  partial  pressures,  is  defined  by 

*  tJLjk  '  U9) 

&  ^  b 
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where  for  any  general  chemical  reaction  between  gases 


&A  +  yy)  M  -f  T)  N+  <  •  . 


(30) 


a,  b,  ...,m,n,  ...  represent  stoichiometric  coefficients 
for  the  species  A,  B,  ...,M,N,  ...  . 

r"  ® 

Kp  is  related  to  the  standard  free  energy  change  Ar  , 

tS*  -  —  RT  K p 

(3D 

where  R  is  the  gas  constant  and  T  is  the  absolute  tempera¬ 
ture. 

The  standard  free  energy  of  formation  of  AI2O3  at  the 

Ignition  temperature,  2300  K,  is  -148  Kcal  per  mole  (Ref  14), 

a  very  negative  value  implying  a  very  stxvng  affinity  of  the 

metal  for  oxygen.  In  fact,  the  reaction  will  proceed  at  a 

-14 

partial  pressure  of  10  atmospheres  of  molecular  oxygen. 

The  kinetics  of  the  chemical  reaction 

4/U  +3  0,  -*  2  AhO, 


are  not  well  known.  But  there  is  reason  to  believe  that  at 
the  temperatures  associated  with  the  diffusion  flame  of 
burning  A1  that  the  reactions 


At  *  O  —  MO 

(33) 

A  AS  + 

(34) 

4-AtO  +  0t—2A<!10j 

(35) 
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all  proceed  (Ref  3*58)  which  suggests  a  chain  reaction  in 
which  the  presence  of  atomic  oxygen  as  a  chain  carrier 
significantly  speeds  the  rate  of  reaction  (Ref  16:366). 
Spectrographic  analyses  of  the  impact  flash  show  the  pres¬ 
sure  of  A10  (Ref  4:33)  reaches  a  maximum  4-180  sec  after 
the  ballistic  shock.  As  a  point  of  interest  Fig  19  indi¬ 
cates  the  partial  pressure  of  atomic  oxygen  as  a  function 

of  velocity.  Note  that  at  8000  ft/sec  the  partial  pressure 

-2 

of  0  in  air  is  10  atmospheres. 
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Backman  (Ref  4:31)  reports  that  by  replacing  air  in 
the  box  with  nitrogen  or  argon  that  the  pressure  rise  Is 
undiminished  when  compared  to  air  in  the  box.  He  concludes, 
logically,  that  burning  cannot  occur.  Assuming  that  he  used 
99£  pure  N2  or  Ar  the  remaining  gas  being  air,  and  approxi¬ 
mately  one-third  of  the  plug  is  shattered  into  small  enough 
pieces,  less  than  100  /x  diameter,  there  is  still  twice  as 
much  02  necessary  to  consume  the  A1  fragments. 

Summary  of  Damage  Mechanism 

In  an  attempt  to-  explain  the  mechanism  for  a  sudden 
rise  in  pressure  in  an  Al  box  when  Impacted  by  A1  fragments 
at  or  above  8000  ft/sec  the  following  hypothesis  was  pre¬ 
sented. 

First,  a  source  of  high  velocity  secondary  fragments 
must  be  found.  This  may  occur  by  a  mechanism  as  postulated 
in  the  Heyda  plug-shattering  model.  For  Al-Al  systems  the 
velocity  where  this  occurs  is  in  the  neighborhood  of  7000- 
12,000  ft/sec  depending  on  geometry. 

Second,  at  8000  ft/sec  the  hypersonic  conditions  pro¬ 
vide  a  source  of  heat  sufficient  to  raise  the  temperature  of 
the  fragments  to  slightly  above  the  ignition  point  of  Al. 

Third,  at  8000  ft/sec  significant  dissociation  of  mole¬ 
cular  oxygen  to  atomic  oxygen  takes  place  which  provides  a 
very  energetic  oxidizer  to  the  system. 

Therefore,  as  every  volunteer  fireman  knows,  if  Ignition 
temperatures,  oxidizer  and  fuel  are  present  a  fire  must  follow. 
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III.  Design  Considerations 

Since  any  application  of  the  synergistic  damage  mech¬ 
anism  would  most  likely  Involve  weapons,  this  section  will 
Indicate  some  factors  which  would  be  considered  In  such  a 
design.  Practically,  the  probability  of  achieving  enhanced 
damage  would  depend  directly  on  creating  a  large  number  of 
micron  size  fragments  moving  in  a  closed,  or  nearly  closed 
box  at  a  velocity  In  the  neighborhood  of  8000  ft/sec  or 
higher.  The  design  of  a  weapon  would  revolve  about  creation 
of  a  large  number  of  suitably  shaped  fragments.  All  designs 
do  not  involve  only  Al-Al  impacts.  A  computer  program  was 
prepared  to  evaluate  the  plug-shattering  phenomenon  for  a 
variety  of  materials. 

Plug-Shattering  for  Other  Than  Al-Al  Impacts 

From  the  Heyda  plug-shattering  model  it  was  concluded 
that  when  the  ratio  of  projectile  radius  to  plate  thickness 
was  greater  than  the  parameter  S,  defined  by  equation  18, 
plug-shattering  may  occur.  Fig  5  and  the  charts  In  Appendix 
B  present  a  relationship  between  the  parameter  S  and  the  ratio 
of  densities  before  and  after  impact  in  the  plate.  Know¬ 
ing  the  density  ratio  it  is  then  possible  to  calculate  the 
particle  velocity  in  the  plate  by  equation  18.  Then  from  the 
material  properties,  specifically  the  emperical  constants  A 
and  B,  a  critical  velocity  Vo  associated  with  incipient  plug- 
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shattering  may  be  determined  from  equation  26.  Note  that  the 
alike  material  impact  is  not  necessary  for  the  plug-shatter¬ 
ing  event  to  occur. 

For  unlike  materials,  for  example,  a  stainless  steel 
projectile  impacting  an  aluminum  plate,  the  '■'litical  velocity 
for  shattering  is  somewhat  less  than  that  of  an  aluminum 
projectile  impacting  a  like  plate.  Fig  20  illustrates  that 
for  a  given  set  of  conditions,  the  stainless  projectile 
causes  shattering  at  considerably  lower  velocity. 


Fig  20.  Velocity  Required  for  A1  and  Steel  Projectiles  at 
Incipient  Plug-Shattering. 
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If  curves  for  a  projectile-plate  system  are  available 
then  the  calculation  of  a  critical  velocity  is  very  straight 
forward.  Appendix  B  presents  curves  plotting  the  relation¬ 
ship  for  a  variety  of  materials  which  may  be  used  for  the 
calculation.  However,  if  curves  for  the  material  of  interest 
are  not  available  the  problem  is  somewhat  more  difficult. 


The  definition  of  the  parameter  S  is  given  by 

s  =  (18) 

By  assuming  a  density  ration,  '#/^  the  particle  velocity 
is  found  immediately  from 


(24) 


and  the  shock  velocity  in  the  plate  D  is  then  found  from 


1):  A+  Bu  <23) 

The  value  for  the  hydrodynamic  sound  speed  C„  is  not 

n 

readily  available  or  calculable  from  elastic  theory  for 
pressure  regimes  of  400  kilobars.  Fortunately,  AFWL-TR-69- 
38,  Compilation  of  Hugoniot  Equations  of  State,  contains  an 
emperical  equation  of  state  for  a  wide  variety  of  materials. 
(It  also  contains  values  for  A  and  B  and  the  density.) 

Recall,  that  the  definition  of  the  sound  speed  is  given 
by 


for  an  isentroplc  wave. 
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The  pressure  In  the  plate  P  associated  with  an  Impact 
velocity  may  be  calculated  by 

(37] 

The  density  ratio  was  used  as  a  starting  point  so 
the  density  of  the  plate  behind  the  shock  wave  Is  simply 


f-  -E 
V  /p. 


(fyp  )«.****,«< 


The  hydrodynamic  sound  speed  Is  then  approximated  by 
taking  a  small  Increment  about  the  pressure  and  calculating 
the  slope  of  the  equation  of  state.  The  assumption  of  an 
lsentroplc  wave  can  be  Justified  by  observing  that  the  sound 
speed  of  Interest  Is  a  dilation  wave  relative  to  its  very 
high  pressure  environment,  behind  the  planar  shock. 


Pressure  Increase  In  a  Closed  System 

A  simple  calculation  of  the  pressure  increase  In  a  box 
Is  necessary  to  gain  insight  into  the  amount  of  A1  which  is 
necessary  to  raise  the  pressure  and  consequently  apply  an 
impulse. 

Consider  a  closed  box  of  one  cubic  foot  volume  into 
which  A1  fragments  are  Introduced  suddenly.  For  this 
analysis  consider  the  box  rigid  and  that  venting  caused  by 
the  perforations  are  negligible.  This  assumption  is  Justified 
the  fact  that  the  fragments  causing  the  perforations  are 
causing  an  air  shock  interaction  at  the  openings  which  would 
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for  the  first  few  milliseconds  of  the  event  essentially  plug 
the  holes.  The  box  contains  initially  air  at  standard  condi¬ 
tions.  The  first  law  of  thermodynamics  for  a  stationary 
system  with  no  external  work  is 

a  =  -^P  Mc  (39) 

where  A is  the  change  of  enthalpy  for  air  between  the 
initial  and  final  states,  Snip  is  the  mass  of  A1  burned  in 
air  with  a  heat  of  combustion,  .  Substituting  the  defi¬ 
nition  of  specific  heat  at  constant  pressure  gives 


(40) 


assuming  a  constant  specific  heat  for  air  and  an  ideal  gas 
at  constant  volume  gives 

A  9  .  in?  Mc  R 

Ar- - l. -  -  (4i) 

Cpfc 

Some  Justification  is  required  for  the  ideal  gas  as¬ 
sumption  used  in  this  analysis.  The  first  reason  of  course 
is  the  fact  that  it  is  the  air  in  the  box  which  is  being 
heated  and  compressed.  The  second  reason  is  that  a  detailed 
model  of  a  similar  situation  was  constructed  assuming  a 
statistical-mechanical  model  of  the  gas  (Ref  3*103)  with  the 
result  that  the  exact  solutions  differed  only  a  few  percent 
from  that  assuming  ideal  gases  and  neglecting  Intermediate 
products. 

Figure  21  represents  the  increase  in  pressure  per  cubic 
foot  of  an  air  filled  box  as  a  function  of  the  amount  of  A1 
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MASS  OF  Al  CONSUMED  ,  mg 


Fig  21.  Pressure  Change  per  Milligram  of  Al  Oxidized 
in  a  One  Cubic  Foot  Volume. 


consumed  by  oxidation. 


Other  Design  Considerations 

The  mechanism  for  sudden  pressure  Increase  Is  a  very 
complex  combination  of  different  phenomena  each  contributing 
to  some  portion  of  the  synergistic  damage  event.  To  analyti 
cally  prescribe  design  criteria  would  be  a  formidable  task, 
at  least.  It  Is  clear  that  the  formation  of  numerous  mlcro- 
scoplc-slze  hlgh-veloclty  secondary  fragments  Is  essential 
and  that  above  8000  ft/sec  for  these  microscopic  fragments, 
combustion  Is  almost  a  certainty.  As  shown  In  Figure  21,  a 
very  small  amou >t  of  Al,  or  other  pyrophoric  metal,  can 
create  significant  pressure  increase  In  a  small  box. 
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It  is  believed  that  great  caution  must  be  observed 
before  attempting  to  apply  this  phenomenon  as  a  damage 
mechanism  for  two  basic  reasons.  First,  if  the  volume  en¬ 
closed  behind  the  leading,  or  punctured  plate,  is  large 
enough  that  numerous  sources  of  secondary  fragments  must  be 
generated,  which  implies  that  a  high  area  density  of  impact¬ 
ing  fragments  is  required.  The  amount  would  vary  with  the 
volume  geometry.  Secondly,  if  the  energy  released  by  the 
combustion  of  secondary  fragments  is  considered  to  be  deposited 
in  very  short  time  at  a  point  and  that  Taylor's  strong  blast 

wave  theori'  applies,  then  the  damage  Impulse  will  be  at- 
-3 

tenuated  as  R  ,  where  R  is  the  distance  from  the  edge  of 
the  energy  source.  Again  there  is  cause  for  great  attenua¬ 
tion  of  a  shock  wave  generated  by  the  impact. 

For  the  above  two  reasons,  and  confirmed  by  experimental 
evidence  for  explosive  oxidations  (Ref  3s 94)  at  much  higher 
velocl;les,  one  would  not  seek  to  protect  the  box  by  insert¬ 
ing  passive  energy  absorbing  material  as  it  will  only  serve 
to  contain  the  free  expansion  of  the  combustion  gases. 
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IV.  Conclusions  and  Recommendations 

Cause  of  Synergistic  Damage  Mechanism 

The  fact  that  there  will  be  an  explosive  oxidation 
following  a  hypervelocity  impact  is  well  known  (Ref  3). 

But  the  Impact  phenomena  occurring  at  intermediate  veloci¬ 
ties  is  very  much  more  complex.  The  cause  of  pret?surb  build¬ 
up  in  a  box  can  be  attributed  to  several  factors,  any  or 
all  occurring  simultaneously.  One  thing  is  reasonably  clear 
however.  The  formation  of  numerous  high-velocity  microscopic- 
sized  secondary  fragments  is  essential  for  any  substantisl 
pressure  increase.  If  the  secondary  particles  are  very 
small  and  have  a  velocity  approximately  8000  ft/sec  then, 
for  aluminum,  at  least,  seme  or  all  will  react  with  the 
available  oxygen  in  the  system  releasing  additional  energy. 

Several  sources  have  reported  significant  increases  in 
the  pressure  beginning  with  an  impact  velocity  of  approxi¬ 
mately  8000  ft/sec.  This  is  because  several  coincidences 
occur  simultaneously  when  a  particle  moves  through  air  at 
that  velocity.  They  ares 

(1)  Stagnation  temperatures  around  the  frag¬ 
ment  are  Just  above  the  melting  point  of 

which 

(2)  allows  reaction  of  atomic  oxygen  with 
the  molten  Al.  Significant  dissociation 
of  0>2  begins  at  8000  ft/sec. 
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(3)  At  approximately  8000  ft/sec  there  la 
sufficient  heat  transfer  for  microscopic 
sized  fragments  to  overcome  any  chemical 
threshold  which  may  exist;  and 

(4)  at  approximately  8000  ft/sec  for  most 
experimental  sized  plates  and  projectiles, 
the  Heyda  plug-shattering  model  predicts 
that  a  source  of  these  high-velocity 
microscopic-sized  fragments  is  probable. 

Application  of  the  Svnergl stlc  Damage  F.echanlsm 

The  combustion  of  high  velocity  secondary  fragments 
following  a  moderately  high  velocity  impact  event  has  the 
capability  to  suddenly  increase  the  pressure  within  the 
box.  The  conditions  which  are  necessary  are  restrictive  in 
that  not  every  impact  event  will  cause  a  secondary  explo¬ 
sion.  Basically  only  small,  flat  projectiles,  which  do  not 
have  ideal  ballistic  characteristics,  will  produce  the 
necessary  secondary  particles  which  are  necessary  for  the 
effect  to  occur.  The  volume  of  the  box  has  an  important 
effect  on  the  amount  of  damage  expected. 

Should  one  want  to  build  a  weapon  relying  on  secondary 
damage  as  a  kill  mechanism  there  would  be,  at  this  level  of 
understanding,  definite  limitations  on  the  confidence  he 
would  have  that  a  kill  would  occur.  Offense  conservative 
criteria,  that  is  the  ability  to  predict  confidently  a 
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level  of  damage,  would  lead  one  co  more  predictable  kill 
phenomena  such  as  penetration  or  blast  overpressure. 

It  Is  recommended  that,  until  a  comprehensive  test 
program  Is  conducted,  the  synergistic  damage  mechanism  be 
treated  as  Just  that,  a  decided  but  unpredictable  enhance¬ 
ment  of  known  damage  mechanisms. 

Future  Activity 

Future  efforts  to  understand  this  damage  mechanism 
should  be  based  on  two  primary  objectives. 

First,  confirmation  through  a  systematic  test  program 
of  the  Heyda  plug-shattering  model  Is  required.  Emperlcal 
data  which  would  confirm  or  deny  the  existence  of  plug- 
shattering  is  scarce  and  what  is  available  was  gathered  for 
other  purposes.  An  extension  of  the  model  to  spherical 
projectiles  was  proposed.  This  extension  must  be  validated 
for  other  structural  materials  than  aluminum. 

As  a  corollary  to  this  work,  further  analytical  work 
on  the  effect  of  spacing  of  multiple  fragments  on  a  plate 
is  necessary  to  develop  additional  mechanisms  for  the  genera¬ 
tion  of  high-velocity  secondary  fragments. 

A  second  area  of  investigation  would  involve  the 
determination  of  the  particle  size  and  velocity  of  the 
secondary  fragments  emerging  from  the  impact  event.  In 
this  study,  only  heuristic  arguments  could  be  presented  on 
the  determination  of  the  size  and  velocity  of  the  particles 
following  the  convergence  of  rarefaction  waves  in  the  plug. 
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This  investigation  would  logically  be  a  part  of  the  verifi¬ 
cation  of  the  Heyda  model  verification. 
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Application  of  Heyda  Plug-Shattering  Model  for  Calculation 
of  Impact  Velocity  for  Incipient  Shattering 


For  test  data  reported  In  AFATL-TR-70-112  (Ref  8 s 146 ( 


Plate  Thickness,  T: 


.635  in,  Steel 


Projectile  radius,  Rs  .391  In,  Steel 
The  critical  condition  for  Incipient  plug- shattering  is 

5  r  *17*  , 

Sc  .(.**/. 31,  =  l.u  4 

From  Fig  5f  for  steel,  when  S=  -  -5o2. 

The  particle  velocity,  U  ,  in  the  plate  is 

u  =  — 

/-%  ~  - 

The  critical  impact  velocity,  Vo  Is 

\  “  _ _ _ 

Vo  ~  u  +  J  A  u  (/}  4  6  u )  -  A  1 

For  like-material  impact  equation  (26)  simplifies  to 


V0~-.  u  +  ~  Ja*  +  4  ABu+  4Z 


•— 1 

^  -  A 

_  1 


(26a) 


From  AFWL-TR-69-38  (Ref  10:79),  values  for  constant  A  and 


£  are: 


A*  Z.m  ,  S  ^  b9 10 


Then  [A  =  •  3  9?  3 
and  from  equation  26a, 


cm^sec,  from  equation  24 


)/ --  cm^asec- ;26/?  ft/sec, 
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Prediction  of  Incipient  Shattering  by  Hey da  Plug-Shattering 
Model  for  a  Spherical  Projectile 


From  test  data  reported  In  Hopkins,  Lee  and  Swift  (Ref  13* 
343) 

Plate  Thickness,  T:  O.787  mm,  606I-T6  A1 

Projectile  diameter,  d:  3.17  ram,  Aluminum, 

From  Fig  8,  for  Incipient  shattering,  let  V0  -  3.0  Km/sec 
Values  for  constants  are  from  Ref  16:25. 

A-  .sott 

6=  /Vzo 


f.. 


^.?3  3  pm/cc 


If  the  projectile  were  cylindrical  the  parameter  S  which  Is 
equal  to  T/e  for  Incipient  shattering  Is 


5  -  V. 8  - 


111 

/.SIS’ 


w?7 


Since  there  Is  no  value  for  S  as  low  as  0.497  and  ex¬ 
perience  Indicates  that  Is  usually  between  0.75  and  0.9 

with  a  corresponding  S  approximately  1.6,  the  approach  will 

R* 

be  to  find  an  equivalent  R*.  If  the  ratio  R  Is  approximate¬ 
ly  constant  for  several  geometries  of  lmpact?then  hypothe¬ 
size  a  corrective  factor. 

The  solution  procedure  is  to  substitute  values  into 
equation  26a  and  solve  for  the  target  particle  velocity 

V."  46*1^  -Aj  (26a) 


54 


GAW/MC/74-10 


3.0  -  U  4  ”7"  I  +  4(\hxoS\a[cm 

A(hj.)  Vv 


which  yields 

U  =  I  Km/sec. 

From  equation  24,  solving  for  ratio  — 

P  ^ 

Al .  I _ n, 

?  1  /i+  Biv 

Substituting  values, 

4l,  , _ - 

(  o'W/  4  I  *2’1) 

?* 

7  =  -77/ 

Referring  to  Pig  5,  the  associated  S  is  I.65. 
For  incipient  shatter 


,JV 

/. 


The  ratio  o* 
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The  following  charts  present  the  values  of  the  para¬ 
meter  S  as  a  function  of  the  density  ratio  following  a 
planar  impact.  When  S  is  equal  to  T/R,  where  T  is  the 
plate  thickness,  R  the  cylindrical  projectile  radius 
incipient  plug-shattering  may  be  expected.  The  range  of 
presented  is  believed,  without  experimental  verifica¬ 
tion,  to  be  the  range  of  practical  application. 

The  charts  may  be  entered  either  through  the  system 
geometry,  S,  and  then  reading  directly  the  density  ratio. 
The  curve  represents  only  the  plate  material.  Below  the 
chart  the  velocity  at  impact  of  a  variety  projectile 
materials  is  presented. 
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